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ABSTRACT 

Recent observations of the Crab pulsar by the Energetic Gamma-Ray Experiment 
Telescope (EGRET) on the Compton Gamma Ray Observatory show that the high- 
energy gamma-ray light curve has changed little over the lifetime of the instrument. 
Previous data collected by SAS-2 and COS-B in the years 1972-82, along with earlier 
EGRET data, suggested a 14 year sinusoidal variation in the flux ratio between the first 
and second peaks. The new data from EGRET indicate that the flux ratio is constant. 
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1. Introduction 

High energy gamma ray emission from the 
Crab pulsar was observed by satellite-borne tele- 
scopes for 15 years: in 1972-73 by SAS-2 (Kniffen 
et al. 1974), from 1975-82 by COS-B (Clear et 
al. 1987), and since 1991 by EGRET (Nolan et 
al. 1993; Ramanamurthy et al. 1995). Early ob- 
servations showed possible sinusoidal variation in 
the relative intensities of the two peaks (Wills et 
al. 1982) with a time scale of ~ 14 years, and it 
was suggested that this variation might be due 
to the precession or free nutation of the neutron 
star (Kanbach 1990, Ozel 1991). An apparent 
confirmation of the sinusoidal signal was seen in 
the low energy gamma ray emission (Ulmer et 
al. 1994), matched in phase and period with 
the high energy results, but with a smaller am- 
plitude. EGRET data from 1991 through early 
1994 were consistent with the expected variation 
(Nolan et al. 1993, Ramanamurthy et al. 1995), 
although these observations spanned a time when 
the ratio of the peaks was predicted to be fairly 
constant, near the minimum of the sinusoid. 

EGRET observations have now extended the 
available data by over two years. The most recent 
data were expected to be 4 — 6a from the average 
of the previous values if the sinusoidal model is 
correct. 

2. Observations and Analysis 

The EGRET instrument is a spark chamber 
gamma-ray telescope with an energy range of 
30 MeV - 30 GeV. Details of the instrument de- 
sign, calibration, and standard analysis software 
are given in Thompson et al. (1993). 

All viewing periods where EGRET was pointed 
within 20° of the Crab were analyzed, with the 
exception of viewing period 0021 (1991 Jul 8- 
15), in which there was a large solar flare. The 
Compton Observatory viewing period numbers 
and dates for these observations are shown in Ta- 
ble 1. The eight viewings numbered 4120 through 
5280 (1995 Feb. through 1996 Aug.) have been 



completed since the time of the previous work of 
Ramanamurthy et al. (1995). 

All photons with measured energy above 50 
MeV, which were within an energy-dependent 
cone of half-angle 6 ma x were used in this analy- 
sis. The angle maxj chosen such that 68% of the 
photons originating from the pulsar are within 
the acceptance cone, is given by (Thompson et 
al. 1993) 

Omax = 5°.85 x (E /lOOMe^ )~ a534 . 

The arrival time of each detected photon was 
transformed to Solar System Barycentric Time 
using the DE200 ephemeris, then binned accord- 
ing to the pulsar phase at that time, determined 
from the Princeton Pulsar Timing Database (Ar- 
zoumanian et al. 1992). This analysis was 
performed using the PULSAR program (Fierro 
1995). 

As seen in Figure 1, the light curve was divided 
into several sections, including Peak 1 (phase .94 
- .04), Peak 2 (phase .32 - .46), and the off-pulse 
Background (phase .46 - .94). These definitions 
follow those used in the COS-B analysis (Wills et 
al. 1982) and are similar to those used by Nolan 
et al. (1993) and Ramanamurthy et al. (1995). 
The background (from the Crab nebula, nearby 
sources, and the diffuse Galactic radiation) was 
assumed constant as a function of pulsar phase. 
The off-pulse count rate was then used to find 
the background-subtracted counts estimates of 
the two peaks (PI and P2). In order to avoid 
any effects of changes in instrument performance, 
the evolution of the ratio P2/P1 was examined 
(as in previous analyses). 

The two peaks have slightly different energy 
spectra (Nolan et al. 1993), and the EGRET re- 
sponse at different energies has changed at differ- 
ent rates (Esposito et al. 1997). Thus the ratio 
P2/P1 is affected by the changes in instrument 
performance over time. Calculations of this ef- 
fect, however, indicate that it is an order of mag- 
nitude smaller than the errors in P2/P1 due to 
Poisson fluctuations. 
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Fig. 1. — Gamma-ray phase histograms of photons with 
E > 50 MeV for EGRET observations of the Crab in 1991- 
2 (a) and 1995-6 (b). The vertical dashed lines indicate 
phase boundaries used in the peak height analysis. The 
horizontal dashed lines indicate the background level as 
determined from the data. 



The differences in the energy responses of 
the SAS-2, COS-B, and EGRET instruments 
are larger than the variation in the EGRET re- 
sponse. However, statistical errors in the previ- 
ous instruments' data are larger as well. Thus 
the value of P2/P1 obtained from EGRET data 
should be comparable to that obtained with SAS- 
2 and COS-B. 

3. Results 

The values of P2/P1 obtained, together with 
the la errors, are shown in Table 1. A reduced 
data set, where nearby points are joined for clar- 
ity, is shown in Figure 2. The data were fit with 
a constant, yielding P2/P1 = .54 ± .03, with 
X 2 = 6.01 with 20 degrees of freedom (DOF). 
Such a low value of x 2 might imply that the er- 
rors in the data were over-estimated. In this case, 
however, the errors arise purely from statistical 
Poisson fluctuations, and the low value must oc- 
cur purely by chance. The data are very consis- 
tent with a constant value of P2/P1. The data 
were also fit with a straight line (x 2 = 5.96 with 
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Fig. 2. — Variation in the ratio of the two peaks in 
the Crab light curve for E > 50 MeV (from Table 1). For 
clarity, the 21 observations are grouped into 6 data points, 
where each point represents the average of several nearby 
observations. Error bars are la. The dashed line is the 
average of all EGRET observations. 



19 DOF), and a quadratic (x 2 = 5.92 with 18 
DOF). Neither result gives a significantly better 
fit: the EGRET data are most consistent with 
no variation. 

The EGRET data were also analyzed in con- 
junction with the SAS-2 (Kanbach 1990) and 
COS-B (Clear et al. 1987) data. The best fit 
sinusoid to the previous instruments' data, 

P2/PI = 0.85-0.56 sin (2tt (T - 1975.67) /13.3) , 

where T is the year of the observation, and to 
the combined data set, 

P2/P1 = 0.544-0.060 sin (2tt (T - 1976.48) /11.55) 

are shown together with the data and the average 
value of P2/P1 in Figure 3. 

As can be seen, the most recent EGRET ob- 
servations (the last data point) are much less con- 
sistent with the large amplitude sinusoid. The 
constant value of 0.528 ± 0.027 gives a x 2 = 24.0 
(with 27 DOF), indicating a good fit. The si- 
nusoid fit to the combined data gives a period 
of 11.6 years with x 2 = 21.2 (with 24 DOF), 
which does not represent a significant improve- 
ment. Thus, using the combined data sets, the 
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Fig. 3. — Variation in the ratio of the two peaks in the 
Crab light curve from SAS-2 (1973), COS-B (1975-1983), 
and EGRET (1991-96), where the EGRET data set has 
been reduced as in Fig 2. The dotted line is the best fit 
sinusoid to the pre-EGRET data. The solid line is the 
best fit sinusoid to all the data. The dashed line is the 
average of all the data. 

data are most consistent with no variation in 
P2/P1. 

The light curves obtained from Phase 1 data 
(Apr 91-Sept 92) and from Phases 4 and 5 (Feb 
95-Aug 96) are shown in Figure 1. The overall 
shape seems to have changed little, in contrast 
with the expected change if the 14 year cycle were 
correct. A x 2 test was performed to compare the 
two light curves in a quantitative way. In order 
to take out possible systematic effects, both an 
additive offset and a multiplicative factor for the 
second light curve were fit to the data. The re- 
sulting x 2 = 58.5 (with 48 DOF) is consistent 
with no change in the light curve. Without the 
offset, x 2 = 87.9 (with 49 DOF) was obtained, 
indicating an inconsistency at the 99.95% confi- 
dence level. The offset required indicates a lower 
background level in the later observations. This 
could be due to a change in the nebular emis- 
sion (de Jager et al. 1996), or could be a result 
of changes in the performance of EGRET. As 
the gas in the spark chamber ages, the sensitive 
area at low energies decreases, which decreases 
the width of the average point spread function. 
This effect might lower the background in the 



4. Summary and Conclusion 

Recent observations with EGRET have pro- 
vided data relevant to the reported variation in 
the P2/P1 ratio of the Crab pulsar. Data from 
SAS-2 and COS-B suggested a sinusoidal varia- 
tion in this ratio. The EGRET data, both taken 
alone and in conjunction with the data from pre- 
vious instruments, are most consistent with a 
constant value of P2/P1. Examination of the 
light curves from early and later observations 
shows no distinct changes in the pulsar's light 
curve. The EGRET data cannot, of course, rule 
out past variability in the P2/P1 ratio. Future 
observations by EGRET and its successors may 
allow a more precise characterization of the long 
term behavior of the Crab's light curve. 
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ment NCC 5-95 (HSC), NASA Grant NAG5- 
1605(SU), and NASA Contract NAS5-96051 (NGC) 
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Table 1 

Crab P2/P1 for Each EGRET Observation 



Viewing Period 


Dates 




P2/P1 


0002 


Apr 22-28 


1991 


0.57 ±0.11 


0003 


Apr 28-May 1 


1991 


0.54 ± 0.15 


0004 


May 1-4 


1991 


0.65 ±0.18 


0005 


May 4-7 


1991 


0.55 ±0.14 


0010 


May 16-30 


1991 


0.52 ±0.07 


0360 


Aug 11-12 


1992 


0.66 ± 0.34 


0365 


Aug 12-20 


1992 


0.37 ±0.15 


0390 


Sep 1-17 


1992 


0.65 ±0.17 


2130 


May 23-29 


1993 


0.45 ±0.19 


2210 


May 13-24 


1993 


0.49 ±0.12 


3100 


Dec 1-13 


1993 


0.55 ±0.16 


3211 


Feb 8-15 


1994 


0.59 ±0.11 


3215 


Feb 15-17 


1994 


0.61 ±0.31 


4120 


Feb 28-Mar 7 


1995 


0.52 ±0.15 


4130 


Mar 7-21 


1995 


0.58 ±0.12 


4200 


May 23-Jun 6 


1995 


0.45 ±0.15 


4260 


Aug 8-22 


1995 


0.50 ±0.20 


5020 


Oct 17-31 


1995 


0.73 ±0.14 


5260 


Jul 30-Aug 13 


1996 


0.50 ±0.12 


5270 


Aug 13-20 


1996 


0.69 ±0.21 


5280 


Aug 20-27 


1996 


0.48 ± 0.22 
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